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ABSTRACT 

We present VR magnification ratios of QSO 0957+561, wliicli are inferred 
from tlie GLITP light curves of Q0957+561A and new frames taken with the 
2.56m Nordic Optical Telescope about 14 months after the GLITP monitoring. 
To extract the fluxes of the two close components, two different photometric tech- 
niques are used: pho2comC and psfphot. From the two photometric approaches 
and a reasonable range for the time delay in the system (415-430 days), we do not 
obtain achromatic optical continuum ratios, but ratios depending on the wave- 
length. Our final global measurements Am^B(Ay) = 0.077 ± 0.023 mag and 
AmABi^R.) — 0.022 ± 0.013 mag (Icr intervals) are in agreement with the Oslo 
group results (using the same telescope in the same seasons, but another photo- 
metric task and only one time delay of about 416 days). These new measurements 
are consistent with differential extinction in the lens galaxy, the Lyman limit sys- 
tem, the damped Ly a system, or the host galaxy of the QSO. The possible values 
for the differential extinction and the ratio of total to selective extinction in the 
V band are reasonable. Moreover, crude probability arguments suggest that the 
ray paths of the two components cross a similar dusty environment, including a 
network of compact dust clouds and compact dust voids. As an alternative (in 
fact, the usual interpretation of the old ratios), we also try to explain the new ra- 
tios as caused by gravitational microlensing in the deflector. From magnification 
maps for each of the gravitationally lensed images, using different fractions of 
the surface mass density represented by the microlenses, as well as different sizes 
and profiles of the l^-band and i?-band sources, several synthetic distributions of 
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[AmAB(Ay),AmAs(AK)] pairs are derived. In some gravitational scenarios, there 
is an apparent disagreement between the observed pair of ratios and the simu- 
lated distributions. However, several microlensing pictures work well. To decide 
between either extinction, or microlensing, or a mixed scenario (extinction + 
microlensing), new observational and interpretation efforts are required. 

Subject headings: dust, extinction — galaxies: ISM — gravitational lensing — 
quasars: individual (QSO 0957+561) 

1. Introduction 

The first lensed quasar QSO 0957+561 (Walsh, Carswell, & Weymann 1979) is an 
important laboratory to study the evolution of multiwavelength magnification ratios in a 
lens system and the nature of the involved physical phenomena. However, although a 
lot of magnification ratios were measured during the last twenty years, there is no a fair 
and complete picture accounting for them. The magnification ratio (in magnitudes) is de- 
fined as the difference between the magnitude of Q0957+561A and the time delay cor- 
rected magnitude of Q0957+561B. At an observed wavelength A and time t, the ratio is 
Am^^(A,t) = m^(A,t) — mB{\,t + At^^), where AIab is the time delay between the two 
components. Different studies have stablishcd that the radio magnification ratio docs not 
depend on the time and it is close to — 0.31 mag (e.g., Conner, Lehar & Burke 1992; Garrett 
et al. 1994). Moreover, the optical emission-lines ratio agrees with the radio magnification 
ratio: AmAB(radio) Am^B (emission — lines) ^ — 0.31 mag (e.g., Schild & Smith 1991). 
These radio/optical results suggest that the macrolens magnification ratio must be of — 0.31 
mag. 

On the other hand, the measurements in several optical filters (mostly containing optical 
continuum light) disagree with the macrolens ratio. The i?-band magnification ratio has 
been more or less constant for about 15 years, but relatively far from the expected ratio for 
no cxtinction/microlensing (Pelt et al. 1998; Oscoz et al. 2002; Ovaldsen et al. 2003a): 
< ArriABi^R, 1986 — 1999) < 0.1. This ratio was intensively derived from modest (~ Im) 
telescopes in different observatories. Another intriguing result was obtained from the Apache 
Point Observatory (APO) data. Kundic et al. (1997) used the 3.5m telescope at the APO, 
and estimated magnification ratios of AmAB{\, 1995) ~ 0.12 mag and AmAsiK, 1995) ~ 
0.21 mag. These last gr ratios were analyzed by Press & Rybicki (1998), who did not obtain 
any fair conclusion. However we feel that the r-band ratio is strongly biased because of the 
methodology to derive it. Prom the gr fight curves, Kundic et al. (1997) simultaneously 
fitted the magnification ratios (magnitude offsets) and the contaminations of the lens galaxy. 
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Nevertheless, as the contamination of the lens galaxy in a given optical band is basically a 
noisy magnitud offset, the Kundic et al.'s technique only works when the magnification ratio 
is noticeably larger than the contamination. Thus, whereas the (7-band ratio may be close 
to the true ratio in that optical band, the r-band ratio could represent the sum of two 
contributions: a part related to the true magnification ratio and an important contribution 
due to the lens galaxy. We note that the intensive i?-band data agree with this reasonable 
hypothesis. 

Once we solved the R — r paradox, it seems still necessary to do a detailed analysis 
of the magnification ratios in the red and blue regions of the optical continuum, study the 
possible evolution of the ratios, and try the interpretation of the results from different per- 
spectives (e.g., extinction and microlensing) . In order to justify the optical continuum ratios, 
gravitational microlensing scenarios are so far the favorite ones. Two facts determine this 
choice. First, radio-signals and optical emission-lines reasonably come from regions much 
larger than the optical continuum sources, so microlensing exclusively would change the flux 
from the optical continuum (compact) sources. Moreover, a dusty environment could lead 
to emission-lines ratios different to ~ — 0.31 mag. Second, there is ambiguity on the chro- 
maticity/achromaticity and evolution/rest of the ratios (using the oldest data, some authors 
claimed the existence of a drop in 1983: AniABi^R, 1983) ~ —0.3 mag, which would repre- 
sent the recuperation of the macrolens magnification ratio and would be difficult to interpret: 
a void in a dust cloud crossing the A image, or an interface between two consecutive dust 
clouds or microlensing events?). This ambiguity permits to assume achromaticity and rest, 
which are fully consistent with sources moving across a roughly constant magnification pat- 
tern (on the scale and path of the sources). However, Press & Rybicki (1998) remarked that 
the excessive constancy of the optical ratios from 1980 through 1995 (and so, the constancy 
of the magnification) is relatively discordant with a microlensing scenario. 

In this paper (Section 2) we present new VR magnification ratios of QSO 0957-1-561. 

The data correspond to the 2000/2001 seasons. In Section 3, the results are compared 
with the predictions of a dusty intervening system. In Section 4, we study the feasibility of 
an interpretation based on a population of microlenses within the deflector (lens galaxy -|- 
cluster). Finally, the main conclusions and future prospects are included in Section 5. 

2. Optical continuum magnification ratios 

We observed QSO 0957+561 from 2000 February 3 to 2000 March 31, as well as in 
2001 April. All observations were made with the 2.56m Nordic Optical Telescope (NOT) 
in the V and R bands. The main period (in 2000) corresponds to the Gravitational Lenses 



-4- 



International Time Project (GLITP) monitoring, which was already reduced, analyzed and 
discussed (UUan et al. 2003). In 2001 April, we obtained frames at two nights. Unfortunately, 
only the frames on April 10 are useful for photometric tasks. We used the Bessel V and R 
passbands, with effective wavelengths of Xy = 0.536 /im and = 0.645 fim, respectively. 
Therefore, as the source quasar is at redshift Zg = 1.41, we deal with the rest-frame ultraviolet 
continuum (0.222-0.268 /xm). 

From the GLITP brightness records of Q0957+561A and the about 14 months delayed 
Q0957+561B fluxes, one can flnd the V-band and i?-band magnification ratios. Thus, first, 
we take the Q0957+561A fiuxcs in Figs. (6-7) of Ullan et al. (2003). These fiuxcs were 
derived from two photometric approaches: pho2com and psfphot, which led to two datasets 
consistent each other. Second, using the pho2com task, the subtraction of the lens galaxy in 
the Q0957-I-561B region is not perfect, and some additional correction must be incorporated 
(Serra-Ricart et al. 1999). UUan et al. presented GLITP correction laws that qualitatively 
agree with contamination trends for the IAC-80 telescope/camera (operated in the Teide 
Observatory by the Instituto de Astrofisica de Canarias), so we use the pho2com, method 
and the GLITP contamination laws to extract clean fiuxes of Q0957+561B on 2001 April 
10 (at the time Ib)- The combined technique {pho2co'm + corr) is called pho2co'mC. Third, 
we obtain the fiuxes of Q0957-I-561B at through the psfphot task. In this photometric 
technique, the way to determine the brightnesses of the two quasar components and the 
lens galaxy is from PSF fitting. We use the code developed and described by McLeod et al. 
(1998), which was kindly provided us. Several details on the psfphot approach were included 
in section 2 of that paper. The F-band (blue filled circles) and i?-band (red filled circles) 
ratios are presented in Figure 1. While the top panel contains the results from pho2comC, the 
bottom panel includes the psfphot results. The Auiab differences are computed in a simple 
way: given an optical filter, a task and a time delay At^^, we compare the Q0957+561A 
data in a 5-days bin centered on ~ ^t^s and the Q0957+561B fiux at t^- Only relatively 
populated bins arc taken into account, i.e., bins with 3 or more data. It is a clear matter 
that the Am^B(Ay, 2000) = ArriABi^v) and ArriABi^R, 2000) = AmABi^R) similar. 
In other words, given a photometric approach and a time delay, the magnification ratios are 
not achromatic, and it is apparent the existence of a correlation between the ratios and the 
corresponding wavelengths, with the higher ratio for the smaller wavelength. 

However, as we only have one frame in both filters {V and R) at Ib, the VR fiuxes 
of Q0957-I-561B are not so well-determined as the free-of-peaks-of-noise and binned fiuxes 
of Q0957-I-561A. Due to this fact, there is a small discrepancy between the VR ratios from 
the two photometries. Fortunately, Ovaldsen et al. (2003b) reported VR fiuxes of QSO 
0957+561 through five nights (one in 2000 January and four in 2001 March) of intensive 
monitoring at the NOT. The Oslo group used the same telescope in the same seasons. 
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although they reduced the frames from another method (different to pho2comC or psfphot) 
and only tested the VR magnification ratios for a delay of about 416 days. Their typical 
values ArriABi^v) = 0.065 mag (blue open circles in Fig. 1) and AmABi^R) = 0.026 mag 
(red open circles in Fig. 1) agree with our results, and roughly represent a gap larger than 
the gap associated with psfphot, but smaller than the pho2comC gap. Once we confirm 
the absence of strongly biased results, as the VR estimates for a given photometry do not 
seem to be very sensitive to the delay value, we introduce delay-averaged VR ratios. In 
this first step, the effective ratios are: Am^B(Ay) = 0.097 ± 0.011 mag and AmABi^R) — 
0.012 ± 0.010 mag, and AmAB(Ay) = 0.057 ± 0.011 mag and AmAB{XR) = 0.032 ± 0.007 
mag, from pho2comC and psfphot, respectively. The effective measurements are obtained 
from the averages of the typical values and errors for the different delays. In a second step, 
taking into account the estimates from the two photometric techniques, we infer the final 
global measurements. These VR measurements are: ArriABi^v) = 0.077 ± 0.023 mag and 
ArriABi^R) = 0.022 ± 0.013 mag {la intervals). Each la interval accounts for both the 
scatter between the typical values and the formal uncertainty of the photometric methods. 
The two individual contributions are added in quadrature. 



3. Dust system 

We are going to discuss the feasibility of two physical scenarios, which are usually 
introduced to explain optical ratios in disagreement with the expected macrolens ratio. The 
relative extinction of Q0957+561A is our first option. This relative extinction may be caused 
by dust in an intervening system. What system?. As the Milky Way is a bad candidate (the 
psfphot measurements and a hypothetical dust system at 2; = lead to unphysical extinction 
parameters), one can discard the differential extinction inside the Galaxy. However, in QSO 
0957+561, we have four more candidates. Lens galaxies can generate differential extinction 
between the components of lensed quasars (e.g., Falco et al. 1999, Motta et al. 2002). Thus, 
the z = 0.36 lens galaxy is an obvious candidate. It was discovered by Stockton (1980), 
and it is a giant elliptical galaxy residing in a cluster of galaxies. The light related to the 
A component and the light associated with the B component cross over two galaxy regions 
that are separated by ~ 20 kpc, where the Hubble constant is given by if = 100 h km 

Mpc~^, Q — 0.3, and A = 0.7. Besides the lens galaxy, there are two Ly absorption-hne 
systems at redshifts z — 1.1249 and z — 1.3911 (e.g., Michalitsianos et al. 1997). The object 
closer to the observer {z — 1.1249) is a Lyman limit system. For this system, spectra of QSO 
0957+561 showed that the absorption in component A is marginally greater compared with 
component B. Therefore, assuming that more gas implies more dust, a relative extinction of 
Q0957+561A could be produced at ^ = 1.1249. Michalitsianos et al. (1997) did a detailed 
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study of the damped Ly a system oX z — 1.3911. In this far object, there is a strong evidence 
for different H I column densities between the lensed components, with the hne of sight to 
image A intersecting a larger column density. The result iV^(H I) > I) was also 

claimed by Zuo et al. (1997), who suggested the possibility of a differential reddening by 
dust grains in the damped Ly a absorber. At 2; = 1.3911, the ray paths of the lensed images 
are separated by ~ 160 pc (the separation between the ray paths is estimated from the 
scheme by Smette et al. 1992, but using an Q = 0.3 and A = 0.7 cosmology). Apart from 
the lens galaxy and the Ly absorption-line systems, one can also consider the source's host 
galaxy at 2; = 1.41. 

Neglecting gravitational microlensing effects and assuming /2y-dependent Cardelli, Clay- 
ton, & Mathis (1989) extinction laws which are identical for both components, the dusty 
scenario leads to 

AmAB(A) - AmAB(oo) + ^E{B - V){a{x)Ry + h{x)\ (1) 

where Am^s(cxD) is the macrolens ratio or the ratio at A = 00 (in our case, iSmAB^P^^ ~ — 
0.31 mag; see Introduction), and AE{B — V) is the differential extinction. Here, Ry — 
Av/E{B — V) is the ratio of total to selective extinction in the V band, and a{x) and 
b{x) are known functions of a; = (1 + Zdust)/^ (Cardelli et al. 1989). We take the typical 
delay-averaged VR magnification ratios from pho2comC and psfphot, which are considered 
as two extreme cases: a big gap of 0.085 mag {pho2comC) and a small gap of 0.025 mag 
{psfphot). The Ovaldsen et al.'s measurement would correspond to an intermediate case of 
Am^s(Av) — AniABi^R) ~ 0.040 mag. For each photometry, we estimate the two extinction 
parameters, AE{B — V) and Ry, as a function of the redshift of an intervening object, Zdust- 
The results are showed in Figure 2. The sohd line represents the extinction parameters from 
the pho2comC ratios, whereas the dashed line traces the parameters from the psfphot ratios. 
In Fig. 2, the red points are the solutions for Zdust — 0.36 (lens galaxy), the green points are 
the solutions for Zdust = 1.1249 (Lyman limit system), the blue points are associated with 
Zdust = 1-3911 (damped Ly a system), and the yellow points are associated with Zdust = 1-41 
(host galaxy). If we focus on the four candidates (lens galaxy, Ly absorption-line systems, 
and QSO's host galaxy), global intervals AE{B — V) = 0.03-0.10 mag and Ry = 2-9 are 
derived. In the Galaxy, typical paths have a ratio of total to selective extinction Ry ~ 3.1, 
while paths through typical and non-typical extinction regions have values in the range 2-6. 
The Ry values in far galaxies can be even higher. For example, the z — 0.68 spiral lens 
B0218+357 has Ry 7-12 (Falco et al. 1999; Munoz et al. 2004). Therefore, the results 
on Ry seem reasonable. On the other hand, the range for differential extinction agrees with 
typical values in lens systems, and as a general conclusion, the observations can be explained 
by means of a dusty scenario. 
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How is the structure of the dust system?. The ray paths of the two components cross 
either different environments (i.e., B crosses an empty zone and A crosses a dusty zone) or 
the same dusty environment. In the first case, as the emission-hnes sources do not experience 
(important) differential extinction and the optical continuum compact sources suffer it, the 
A zone would include a compact dust cloud in the trajectory of the A beam together with a 
random distribution of similar clumps of dust. The distribution of clouds cannot significantly 
perturb the radiation from the emission-lines sources. If the projected (into the source plane) 
cloud radius is Rdoud — 0.01 pc, i.e., similar to the expected radial size of the whole accretion 
disk around a 10^ M© black hole, then the time of cloud crossing by the optical continuum 
region (OCR) is of about 30 years (considering Rqcr < Rdoud and an effective transverse 
velocity of 600 km s~^; e.g., Wambsganss et al. 2000). If we also consider the existence of A'^ 
dust clouds in front of the broad-line region (BLR) as well as a typical radius Rblr = 0.1 
pc, the surface covering factor of the dust would be f dust/BLR = N {Rdoud / Rblr)'^ = 10~'^N. 
As f dust /BLR should be much smaller than 1 {fdust/BLR « 1), a constraint on the surface 
density of clouds can be derived in a simple way. It must be much smaller than IO^/tt clouds 
pc~^, so the probabihty to see a cloud crossing phenomenon is equal to or less than 10%. 
The situation is better in the second case, when the A and B zones are within the same 
environment. While a homogeneous distribution of dust is not consistent with the optical 
continuum ratios, an inhomogeneous distribution works well. We assume an environment 
constituted by a network of compact dust clouds and compact empty regions {Rdoud — Rvoid 
= 0.01 pc), with a cloud in the trajectory of the A beam and a void in the trajectory of 
the other (B) beam. In this picture, the optical continuum ratios would be anomalous, 
but the emission-line ratios would be very close to the macrolens ratio, because two large 
regions contain a very similar amount of dust. Moreover, both the probability of a cloud 
crossing (component A) and the probability of a void crossing (component B) are of about 
50%. Hence, the probability of a cloud(A)/void(B) crossing is crudely of 25%. For the farest 
candidate (host galaxy), the common environment is a very reasonable approach, and for 
the nearest candidate (lens galaxy), the common environment could be due to the presence 
of an inhomogeneous dust lane with a length of > 20 kpc. 



4. Microlensing 

To explain the chromatic ratios, the gravitational microlensing of QSO 0957+561 is 
our second option. While the emission-lines regions are not magnified from gravitational 
microlensing, the optical continuum compact sources could be affected by this physical phe- 
nomenon. Using the values for the normalized surface mass density (convergence) and the 
shear, we produce 2-dimensional magnification maps with a ray shooting technique (Wamb- 
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sganss 1990, 1999). These values are: ka = 0.22, 7^ = 0.17, kb = 1.24, and 75 = 0.9 (Pelt 
et al. 1998). With respect to the fraction of the surface mass density represented by the 
microlcnses, we consider three cases: (a) all the mass is assumed to be in granular form 
(compact objects), (b) 50% of the mass in granular form, and (c) the mass is dominated 
by a smoothly distributed material, with only 25% of the mass in compact objects. All the 
microlensing objects are assumed to have a similar mass M and are distributed randomly 
over the lens plane. We obtain detailed results for M — IMq, and comment the expected 
results for smaller masses (M < IMq). The effect of the source is taken into account by 
convolving a selected source intensity profile with the magnification maps. We consider two 
different circularly-symmetric profiles: a p = 3/2 power-law (PL) profile and a Gaussian 
(GS) profile. The PL profile was introduced by Shalyapin (2001), and it is roughly consistent 
with the standard accretion disk. The GS profile is not so smooth as the PL one, and it 
was used in almost all previous microlensing studies in QSO 0957+561. In our model we 
assume that the V-band source has a characteristic length (intensity distribution) of Ry = 
3 X 10^"^ cm (small source), Ry = 10^^ cm (intermediate source), or Ry = 3 x 10^^ cm 
(large source), and the source size ratio q — Ry/Rn is given by about 0.8 or 1/3. For the 
standard accretion disk model, one has the relationship R{X) oc A^/^ (e.g., Shalyapin et al. 
2002), and consequently, the standard ratio q ^ 0.8, which is close to the critical ratio (g 
= 1). It is also assumed that the source, the defiector and the observer arc embedded in a 
standard flat cosmology: H ^ 100 h km s"^ Mpc-^ h = 2/3, Q = 0.3, A = 0.7. 

In other recent analyses about QSO 0957+561 microlensing (Refsdal et al. 2000; Wamb- 
sganss et al. 2000), the authors concentrated on the time evolution of Attiab in a given optical 
band. Refsdal et al. (2000) used a GS source and found 90% confidence level constraints 
for the microlenses mass M and the source radius R: M = 2 x 10^^ - 0.5 Mq, R < 6 x 
10^^ cm. On the other hand, Wambsganss et al. (2000) claimed that the microlenses mass 
must be > 0.1 Mq for small GS sources (3o" limit). All these previous constraints arc based 
on the hypothesis of a Gaussian source, and moreover, as the analyses were focused on the 
time domain, some hypotheses about both the effective transverse velocity and the random 
microlens motion were needed. However, our approach is different, because we wish to study 
the spectral evolution of Attiab at a given emission time. In this paper, we discuss both 
the GS and PL profiles and do not need to introduce a dynamical hypothesis. In fact, in 
relation to multiepoch simulations, multiband simulations have the advantage of being free 
from assumptions on the peculiar motions of the quasar and defiector, and the orbits of the 
microlenses in the lens galaxy and cluster centre. 

The key idea is that the VR magnification ratios are generated by two circular concentric 
sources with different size. The common center of the y-band and i?-band sources is placed 
at two arbitrary pixels of the two magnification patterns, and we look for the corresponding 
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[AmAB(Ay),AmAs(Ai?)] pair. For a given set of physical parameters, we can test a very large 
number of pairs of pixels (10^), and thus, obtain a distribution of pairs of magnification ratios. 
In the [AmAB(Av),AmyiB(A_R)] plane, for g = 1 (critical source size ratio), the simulated 
distributions would lie exactly on the ArriABi^R) = Am^B(Av') straight line (critical line). 
However, when q ^ 1, we must infer 2D distributions, including points at non-zero distances 
(in magnitudes) from the critical line. The 2D distributions in the six panels of Fig. 3 have 
been derived by using a deflector with only granular matter (M = IMq) and GS sources. In 
Fig. 3 (and also Fig.4), the results for the small, intermediate and large F-band sources are 
depicted in the top, middle and bottom panels, respectively. The left panels correspond to 
a standard q, and the right panels correspond to a non-standard q {q — 1/3). We also draw 
the rectangle associated with the la measurements. It is apparent that larger differences 
between Ry and imply wider distributions. In Fig. 3, we only see a very interesting 
region of the whole plane, which includes the observed magnification ratios. The reader can 
extrapolate the six behaviours in this region towards a larger region along and around the 
critical line. In order to get a good resolution, instead of the whole 2D distributions, we 
decided in favour of a smaller (but representative) region. What about the PL intensity 
profile?. In Fig. 4, we again consider all the mass in compact objects of M = IMq, but PL 
sources. The new distributions are different to their partners in Fig. 3 (with similar values 
of Rv and Rr). However, for obtaining quantitative conclusions, some additional analyses 
are needful. 

We make 12 histograms N{d) for the 12 scenarios where 100% of the mass is contained 
in compact objects (M = IMq). Each histogram represents the number of points at different 
distances d (in magnitudes) from the critical line. Given the [AmAB(Ay),AmyiB(A_R)] pair, 
the distance from the critical line, d = [Am^B(Ay) — ArriABi^R)]/ V^, ^^e difference 
between ratios, D = AmABi^v) — AmABi^R), differ only in a factor \/2. Therefore, we 
analyze histograms N{d), which are basically D distributions. The points above the critical 
line have negative distances {d < 0), whereas the points below the critical line have positive 
distances {d > 0). All the histograms are showed in Fig. 5: GS profile and standard q (top 
left panel), GS profile and q — 1/3 (top right panel), PL profile and standard q (bottom 
left panel), and PL profile and g = 1/3 (bottom right panel). We use blue lines for the 
small K-band source, green lines for the intermediate K-band source, and red lines for the 
large V^-band source. The four arrows indicate the distances from the critical line to the 
four vertexes of the la rectangle in Figs. 3 and 4. If the true [Am^B(Ay),Am^B(A_R)] pair 
would be within the la box, then the true distance would verify d > 0.013 mag. Here, we 
assume that the true pair is placed at d > 0.013 mag, so probabilities P{d < 0.013 mag) are 
computed. Some probability distributions (histograms) lead to very high values of P{d < 
0.013 mag), and thus, we can rule out the corresponding scenarios. On the contrary, other 
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N — d relationships are relatively consistent with the constraint. The estimates of P{d < 
0.013 mag) appear in Table 1, and we see several values < 90%, which suggest realistic 
physical parameters, or in other words, viable scenarios. In the GS profile case, the feasible 
pictures arc related to either the large \/-band source or an intermediate l^-band source 
together with a non-standard i?-band companion. In the PL profile case, from the P < 
90% criterion, all the scenarios seem to be plausible. 

In Fig. 6, we see the magnification maps for the Q0957+561A (left panel) and Q0957+561B 
(right panel) components as well as the 10000 pairs of pixels (white marks) used to produce 
Figs. 3, 4 and 5. For the B component, the granular mass density is high and the caustic 
network almost takes up the whole map. The 10^ highlighted pixels in each map are ran- 
domly distributed, so the sampling is not biased. In Fig. 7, the chromatic microlensing is 
apparent. The caustic maps in Fig. 6 are convolved with the largest Gaussian sources {Ry — 
3 X 10^^ cm, Rr = 3 X Ry). In the three top panels of Fig. 7, we include the ^-convolved 
map for the A component (left), the i?-convolvcd map for the A component (centre), and 
the difference between the two convolutions (right). The results for the B component appear 
in the three bottom panels of Fig. 7: V^-convolution (left), i?-convolution (centre), and 
difference (right). It is clear that the convolutions with the V-band source are different to 
the convolutions with the i?-band source. This fact gives rise to the appearance of chromatic 
effects, i.e., there are difference signals in the two right panels. 

When all the mass is not contained in compact objects, i.e., a certain mass fraction 
is smoothly distributed, we can repeat the above-mentioned procedure and obtain new 
2D distributions, histograms, probabilities, etc. To avoid an "excess of complementary 
material" (figures), only the P{d < 0.013 mag) values are included in Table 1. Two different 
cases are analyzed in this paper: 50% of the mass in granular form (50% of the mass in 
a smoothly distributed component), and 25% of the mass in compact objects (75% of the 
mass is contained in a smoothly distributed material) . As the percentage of mass in compact 
objects decreases, it is more difficult to find scenarios in agreement with the observational 
constraint. In fact, when the mass is dominated by a smoothly distributed material, almost 
all the physical pictures are ruled out by the observations at the 90% level. Only two rare 
scenarios (where the source sizes are in a relation 1:3) lead to P{d < 0.013 mag) < 90%, 
or equivalently, P{d > 0.013 mag) > 10%. Although from a population of microlenses of 
M — IMq we are able to test an important set of physical situations, we must also think 
about the expected results from a population of microlenses of M < IMq. In the extreme 
(but still plausible) case of M = IO^^Mq, the physical size of the magnification maps (the 
maps are 2048 pixels a side which cover a physical size of 16 Einstein radii) reduces in a 
factor 10 (just the factor that relates the Einstein radii for 1 Mq and 1O~^M0). Now, the 
l^-band and i?-band sources with radii Rr and Ry work as the sources with radii 10 xRr 



- 11 - 



and 10 xRv in the old case (M = IM©). For example, the results in Table 1 for Ry = 3 x 
10^^ cm would be the results for = 3 x lO^'' cm and M = IQ-^Mq. 

We remark that the discussion is based on arrays of 2048 by 2048 pixels, which cover 16 
Re X 16 Re regions. Similar arrays were used by Schmidt & Wambsganss (1998) in the cases 
M > 1O~^M0, and we think that our results are rehable. However, in very detailed studies to 
discriminate between different microlensing pictures, magnification maps with larger spatial 
coverage and better spatial resolution may be useful. Some of these improvements are related 
to the availability of high performance computers. 



5. Conclusions and future work 

We present in this work a very robust estimation of the VR magnification ratios of 

QSO 0957+561. The new measurements are based on observations (in 2000-2001) with 
the 2.56m Nordic Optical Telescope, two different photometric techniques, and a reasonable 
interval for the time delay in the system. Our la estimates ArriABi^v) = 0.077 ± 0.023 
mag and Amj^Bi^R) = 0.022 ± 0.013 mag are supported by the results in Ovaldsen et al. 
(2003b). A first conclusion is the relation Am^B(Ay) > AmABi^R), which means that the 
optical continuum magnification ratios are not achromatic, and moreover, the higher ratio 
corresponds to the smaller wavelength. 

To explain the measurements, two alternatives are explored: a dust system between the 
quasar and the observer, or a population of microlenses in the defiector (the most popular 
picture). We find that the observed ratios are consistent with both alternatives, i.e., differ- 
ential extinction in a dust system and gravitational microlensing in the defiector. We cannot 
rule out the existence of a dusty scenario constituted by a network of compact dust clouds 
and compact empty regions, with a cloud in the line of sight to the Q0957+561A component 
and a void in the trajectory of the Q0957+561B beam. The new estimates cover a narrow 
spectral range (0.536-0.645 /um), and this fact plays a dramatic role in the discrimination of 
possible physical scenarios. We think that future accurate measurements at a large collection 
of wavelengths, with a good coverage of the 0.2-1 /xm range, will decide on the feasibility 
of a dust system in QSO 0957+561 and other things. For example, the GTC (a segmented 
10.4 meter telescope to be installed in one of the best sites of the Northern Hemisphere: 
the Roque de los Muchachos Observatory, La Palma, Canary Islands, Spain. First fight 
is planed for 2005) will have an instrument (OSIRIS) with very interesting preformances. 
OSIRIS will incorporate tunable filters with FWHM from 10^*^ to 7 x 10~^ //m over the 
whole optical wavelength range, and thus, tunable imaging of the two quasar components 
at two epochs (separated by the time delay) could lead to solve the current puzzle (dust or 
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microlenses?) and infer tight constraints on the physical parameters of the favoured alter- 
native. Observations in the ultraviolet wavelength range (< 0.36yum) may also be a decisive 
tool. In particular, if the actual origin of the magnification ratio anomalies is a dust system 
in the lens galaxy 

If the 0.2175 fxm extinction bump would be unambiguously detected (assuming dusty 
scenarios in the 0.36-1.41 redshift interval, the rest-frame wavelength of the peak corresponds 
to 0.30-0.52 yum in the observer frame), then the current extinction/mi crolensing degenera- 
tion would be broken. In this last case (detection of extinction), one may try to derive a very 
rich information: the differential extinction, the ratio of total to selective extinction in the 
V band, the macrolcns ratio, the redshift of the dust, and so on (e.g., Falco et al. 1999). On 
the other hand, if the dusty scenarios would be ruled out from future observations at lots of 
wavelengths, the multiwavelength ratios may be used to decide on the viable gravitational 
microlensing scenarios. Prom a simple analysis, we show that several microlensing scenarios 
are ruled out by the new observations at the 90% level (see Table 1 and Sect. 4), so future 
detailed data could lead to important constraints. The future observational work will require 
new ingredients in the microlensing scheme, e.g., new source intensity profiles (in particular, 
the standard accretion disk profile; Shalyapin et al. 2002), or a two mass population of 
microlenses, describing the granular matter in the lens galaxy and the granular matter in 
the cluster centre (e.g., Schechter, Wambsganss & Lewis 2004). Sophisticated source models 
could be also included. The double ring structure suggested by Schild & Vakuhk (2003) is 
an obvious candidate. These authors have described the i?-band source, but they did not 
quote the chromatic behaviour of the five model parameters: the radius of the large ring, 
the radial thickness of the large ring, the radius of the small ring, the thickness of the small 
ring, and the brightness ratio of the two rings, which is a basic issue to simulate multiwave- 
length magnification ratios. The extinction and microlensing hypotheses can be tested in 
both the frequency and time domains, and several previous efforts focused on the possible 
microlensing variability (e.g., Refsdal et al. 2000; Wambsganss et al. 2000). Obviously, the 
observations in the time domain (evolution of the magnification ratios) have a great interest. 
However, when they are compared to microlensing simulations, the microlensing framework 
must incorporate dynamical information, i.e., the confrontation between the observed evo- 
lution and the microlensing scheme requires additional free parameters (which are needless 
in multiwavelength studies). Finally, we note that the two alternative schemes (only dust 
or only microlenses) might be biased models of the actual situation, and a mixed scenario 
(extinction + microlensing) is also possible. Very recently, to explain spectrophotometric 
observations of the lens system HE 0512—3329, Wucknitz et al. (2003) proposed this mixed 
model. 
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Tabic 1: Probabilities for measuring a pair of magnification ratios (one in the V band and 
tlie other in the R band) at a distance < 0.013 magnitudes from the critical line AmABi^R) 
= ArriABiM)- 



Mass in granular form* (%) 


Intensity profile 


Rv (cm) 


q = Rv/Rr 


P{d < 0.013 mag) (%) 


P < 90%? 


100 


Gaussian 


3 X 1014 


0.8 


93.5 


- 








1/3 


91.4 


- 






IQiB 


fs 0.8 


91.1 


- 








1/3 


81.7 


yes 






3 X 10i5 


« 0.8 


80.0 


yes 








1/3 


68.6 


yes 




Power-law (p = 3/2) 


3 X 10" 


fs 0.8 


89.1 


yes 








1/3 


85.5 


yes 






IQiB 


Ri 0.8 


86.6 


yes 








1/3 


75.3 


yes 






3 X IQis 


^ 0.8 


78.2 


yes 








1/3 


69.0 


yes 


50 


Gaussian 


3 X 1014 


« 0.8 


97.9 


- 








1/3 


96.5 


- 






IQiB 


Ri 0.8 


96.8 


- 








1/3 


90.8 








3 X 10i5 


« 0.8 


91.3 










1/3 


80.4 


yes 




Power-law (p = 3/2) 


3 X 10" 


fs 0.8 


95.6 










1/3 


93.1 








IQiB 


Ri 0.8 


93.9 










1/3 


84.7 


yes 






3 X lO^s 


^ 0.8 


85.9 


yes 








1/3 


75.8 


yes 


25 


Gaussian 


3 X 10" 


0.8 


99.3 










1/3 


99.0 








IQiB 


Ri 0.8 


98.8 










1/3 


96.6 








3 X 10^5 


RJ 0.8 


95.8 










1/3 


89.4 


yes 




Power-law (j> = 3/2) 


3 X 10" 


Ri 0.8 


98.2 










1/3 


97.2 








IQiB 


RJ 0.8 


97.7 










1/3 


92.5 








3 X 1015 


Ri 0.8 


92.3 










1/3 


80.8 


yes 



"All the microlenses are assumed to have a similar mass of M = IMq 
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Fig. 1.— QSO 0957+561 magnification ratios in 2000 (2000/2001 seasons). While the 
top panel contains the results from the photometric method pho2comC, the bottom panel 
includes the results from an alternative photometric task [psfphot). Blue and red points 
are associated with ratios in the V and R bands, respectively. The filled circles are our 
measurements and the open circles are the typical estimates by the Oslo group (they did not 
use pho2comC OT psfphot, but a different task). We are able to find ratios for time delays 
different to 415-416 days. In particular, the 420-430 days interval is very well tested. 
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Fig. 2. — /S.E{B — V) and Ry extinction parameters for dusty objects at 0.36 < Zdust < 1-41. 
The solid line represents the parameters from the typical delay-averaged pho2comC ratios, 
whereas the dashed line traces the parameters from the typical delay-averaged psfphot ratios. 
The red points are the solutions for a dust system within the lens galaxy {zdust = 0.36), the 
green points represent the solutions for Zdust = 1.1249 (Lyman limit system), the blue points 
are related to a dusty damped object {zdust = 1.3911), and the yellow points are associated 
with a dusty host galaxy {zdust = 1-41). 
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Fig. 3. — Synthetic distributions of [AmAs(Ay),AmyiB(Aij)] pairs for a deflector with only 
granular matter (M = IMq) and Gaussian sources. The radii of the K-band sources are Ry 
= 3 X 10^^ cm (top panels), Ry = 10^^ cm (middle panels), and Ry = 3 x 10^^ cm (bottom 
panels). In the left panels, we consider a standard q = Ry/R^ of about 0.8, and in the right 
panels, we take a non-standard source size ratio of g = 1/3. The la measurements in this 
paper are also depicted in each panel (small rectangle). A few simulated pairs are within 
the 1(7 rectangle, so they are fully consistent with the observations. However, we generate 
10^ pairs for each physical picture, with the vast majority of points out of the observational 
box, and even out of the —0.15 < AniABi^v) (in mag) < 0.15 and —0.15 < AniABi^R) (i^ 
mag) < 0.15 region. 



-19- 





Am^giXy) (mag) Am^g{Xy) (mag) 



OD O 



a: o 



J 1 1 1 1 1 1 . 




S.STV* I *i 1 1 1 " 


- 1 1 1 1 1 1 l- 



-0.1 0.1 
Am^jjCX^) (mag) 



QO O 



^o; O 



,^5 



o 

I 



•J •. •!" 



8-1 ■ i .THH-'. ^-^ 



-0.1 0.1 
Am,3(X^) (mag) 



Fig. 4. — Synthetic distributions of [Am^B(Ay),Am^B(Ai?)] pairs for a deflector with only 
granular matter (M = IMq) and p = 3/2 power-law sources. See the caption in Fig. 3 for 
comments on the six panels. 
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Fig. 5. — Histograms N[d) for the scenarios with 100% of the mass in compact objects 
(M = IMq). For a given set of physical parameters, we infer a distribution of 10^ points 
in the [AmAB{^v),^^AB{^R)] plane. In a second step, we obtain the number of points at 
different distances d (in magnitudes) from the critical line AniABi^R) = AmABi^v)- The 
distribution N{d) is computed for points above the critical line {d < 0) and below the critical 
line {d > 0). The results from Gaussian sources and a standard source size ratio (~ 0.8), 
Gaussian sources and a non-standard g(=l/3),p = 3/2 power-law sources and a standard 
q, and p = 3/2 power-law sources and a non-standard q appear in the top left panel, the 
top right panel, the bottom left panel, and the bottom right panel, respectively. We use blue 
lines for Ry = 3 x 10^^ cm, green lines for Ry = 10^^ cm, and red lines for Ry = 3 x 10^^ 
cm. The four arrows indicate the distances from the critical line to the four vertexes of the 
la rectangle in Figs. 3 and 4. 
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Fig. 6. — Details on the caustic maps and the locations of the 10^ pairs of pixels used to 
produce Figs. 3-5. The selected pairs of pixels (white marks) are randomly distributed 
within the magnification maps for the A (left panel) and B (right panel) components. 



Fig. 7. — Example of chromatic microlcnsing. We use the magnification maps in Fig. 6 and 
the largest GS sources {Ry = 3 x 10^'^ cm, Rr = 3 x Ry)- The convolved maps for the 
A component appear in the top panels: ^-convolution (left), i?-convolution (centre), and 
difference between both convolutions (right). On the other hand, the convolved maps for 
the B component are depicted in the bottom panels: ^-convolution (left), i?-convolution 
(centre), and difference (right). The difference signals have a structure that is correlated to 
the structure of the caustic networks. 
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